Collapse and Revival of Entanglement between Qubits Interacting via a Quantum Bus 
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We study the dynamics of the Jaynes-Cummings Model for two level systems (or qubits) interact- 
ing with a quantized single mode electromagnetic cavity (or 'quantum bus'). We show that there is a 
time in between the collapse and revival of Rabi oscillations when the state of the qubit sub-system, 
W) attractor' IS'i'Sely independent of its initial state. This generalizes to many qubits the discovery 
by Gea-Banacloche for the one qubit case. The qubits in such 'attractor' states are not entangled 
either with the field or among themselves, even if they were in the initial state. Subsequently the 
entanglement between the qubits revives. Finally, it is argued that the collapse and revival of en- 
tanglement and the persistence of 'non-classicality' is a generic feature of multiple qubits interacting 
via a 'quantum bus'. 



At the heart of Quantum Information Science there 
is the essential quantum mechanical notion of entangle- 
ment. Measures and dynamics of entanglement are at the 
centre of much current research as there are many fun- 
damental questions yet unanswered [l], Q • Under such 
circumstances the study of simple models which feature 
interesting time evolution of entanglement 1^ , and yet 
are tractable, is of particular importance. In this letter 
we shall highlight one of these. 

Quantum dynamics of two level systems (also known 
as qubits), such as spins in a magnetic field, Rydberg 
atoms or Cooper Pair Boxes, coupled to a single mode of 
an electromagnetic cavity are of considerable interest in 
connection with NMR studies of atomic nuclei , Cavity 

guantum Electrodynamics Q and Quantum Computing 
] respectively. The simplest model which captures the 
salient features of the relevant physics in these fields is 
the Jaynes-Cummings model (JCM) 0] for the one qubit 
case and its generalization for multi qubit systems by 
Tavis and Cummings Q. In this letter we wish to focus 
on the interesting dynamics of entanglement between the 
qubits as described by these models. 

One of the most interesting and surprising predictions 
of the JCM is the 'collapse and revival' of Rabi oscil- 
lations of the occupation probabilities for various qubit 
states as the system evolves, from an initial state which 
is a product of a coherent state |a), for the radiation 
field, and a generic qubit state \ipNq) It is central 

to our present concern that such remarkable dynamics 
occurs only because both the matter and the cavity field 
are treated fully quantum mechanically. Thus, in the 
language of quantum computation we may regard 
the above system as a collection of qubits interacting 
via a quantum bus. Indeed our aim here is to study 
the 'collapse and revival' of entanglement between non- 
interacting qubits induced by a quantum bus 

For clarity let us recall the multi qubit JCM Hamil- 
tonian, for which each qubit labelled i can be either in 
its ground state \gi), with energy Cg^i, or its excited state 



|ei), with energy e^^i. Up to a constant the Hamiltonian 
may be written in the following conventional form 



1=1 1=1 
= - |5»)(5»I, cr+ ^ |e,)(g,|, &r ^ |5,)(e4l) 

where and a are the creation and annihilation opera- 
tors of photons with frequency uj, Mli ~ e^^i — eg,i and Xi 
is the cavity-qubiti coupling constant. Here we consider 
only the case of resonance between the qubits and the 
cavity e.g u = Cli for all i and uniform coupling Xt = X. 

The celebrated 'collapse and revival' can be observed 
in the one qubit case. We define the initial state: 



|*i(0)) = |V^i)|a) 



(2) 



a) = eH"l a = V^e-« and 



where 

iV'i) = (Cglff) + Ce|e)). n is the average number of pho- 
tons in the field. The Rabi oscillations of the probability 
that the qubit is in the initial state at first collapse, on 
a time scale of tc — and then revive at ~ 
This is illustrated in Fig. [T]for Cg = 1, = by plot- 
ting J2'^=o l(5''^l^i(0)l^ where {g,n\ is the state for the 
qubit in the ground state with n photons in the cavity. 

A second notable feature of this time evolution, dis- 
covered by Gea-Banacloche is that at ^tj., \'^i{^tr)) 

again factorises into a qubit part {ip) attractor ^^"^ ^ cavity 
part Moreover, remarkably, the former is given 

by 



/ attractor 



V2 



(e-«|e)±z|5)) 



(3) 



where 9 is the phase of the initial coherent state, for 
all initial conditions such that \Cg\^ + |Ce|^ = 1. Note 
also \'4')attractor attained at t = 3tr/2. Because of 
this strikingly non-linear behavior, following Phoenix and 
Knight [13|], we shall refer to these states, Eq. ([3]), as 
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FIG. 1: (color online) Time evolution for a system with one 
qubit. (a) the entropy of the qubit. (b) the probability of 
being in the qubit's initial state \g). (c) the probability of 
being in the state Wattr actor- the probability of being 

in the 'attractor' state goes to one while the entropy goes to 
zero. The qubit starts in the initial state \g) and the value of 
n = 50. 




t/tr 

FIG. 2: (color online) Time evolution for a system with two 
qubits. (a) the entropy of the qubits. (b) the probability of 
the two qubit state \gg). (c) the probability of being in the 
two qubit 'attractor' state {''P'z)'^^^,.^^^^. when the initial phase 
of the radiation field is ^ = 0. The two qubit 'attractor' 
state is reached at tr/4. The initial state of the qubits is 
-7=(|ee) + \gg)) and the value of n = 50. 



'attractors'. The probability that the qubit is in the state 

\'^^ attractor^ 

as depicted by J^'^^o \{'^tttractor,n\'^ li*))] 
is also shown in Fig. [1] together with the von Neumann 
entropy S"3(t) = -Tr {p'-^{t) lnp'3(t)) associated with the 
reduced density matrix pQ(t) = Tr^ (|^'i(i)> of 
the qubit by tracing over the field. Clearly, at t = ^t,. 
the entropy S'^{t) tends to zero and attains it as n — *■ oo, 
indicating that the radiation field and the qubit are not 
entangled [l^ . 

Prompted by these results we have investigated the 
A^g > 1 qubit evolution, starting in the state l^'w, (0)) ~ 
|'(/)jv,)|Q!), and found that the spin coherent states [ij] 

can also be regarded as 'attractors' in a similar, dynam- 
ical, sense as outlined above. The only difference is that 
in the Nn > 1 case )^^^ ^ will occur only for 

^ 1 ' attractor ^ 

a restricted range of initial conditions which we shall 
call basin of attraction. At the attractor time there is 
no entanglement between the qubits and the radiation 
field, furthermore because iV'Ar,)^^^^^^^^^^ is a- product of 
one qubit states, the qubits are not entangled with each 
other. Below we explore the implications of this ob- 
servation for the dynamics of entanglement between the 
qubits. 

From the above point of view the simplest case of in- 
terest is that of two qubits e.g. Ng = 2. In this case the 
time evolution described by |5'2(<)) is readily found [lit . 
For the most general, normalized, initial state 

IV'2) = Cee |ee) + Ceg \eg) + Cge \ge) + Cgg \gg) (5) 

of the qubit sector the exact analytical solution will be 
given elsewhere [lij . Here we consider only the sector 



determined by the restrictions: a = e*^Cee = e ^^Cgg 

and \J \ ~ = Ceg = Cge- As will be illustrated 
presently these define the 'basin of attraction' for the 
'attractor' \i^2)^ttractor- Namely, for any values of a sat- 
isfying < |a| < 1/3 in 



1^2) - a {e-^0\ee) + e^'\gg)) + ^- - \a\^ {\eg) + \ge)) 

, , (6) 

the probability that the two qubits are in the 

state |V'2)attractor- (gi^en by P2attractor{t) 
{'^iattractor\P^iiMt attractor))^ wiU rCach 1 at SOme 

time t* . An example of such behavior (for = 0) is 
shown in Fig. [S] To highlight the similarity with the 
analogous phenomena in the one qubit case (Fig. [1]) we 
also show the entropy S^{t). This is calculated from 
p*^, the two qubit density matrix reduced with respect 
to the cavity field coordinate, which describes a mixed 
state for most times t. Notably, at t* = where 
P2 attractor [t) = 1, the cutropy S^{t) tends to zero in the 
large n limit, indicating that the system of two qubits is 
not entangled with the field. 

The interesting new feature of the two qubit case as op- 
posed to the one qubit case is that the former is in general 
host to entanglement between qubits and this provides an 
opportunity to study the dynamics of such entanglement. 
For example, whilst almost all of the initial states in the 
'basin of attraction' given in Eq. ([6]) describe entangled 
qubits they all evolve into \'^2)~^ttractor at i = \tr where 
they are not entangled. We plot the pure state tangle of 
the initial condition defined as r = 4 \ CeeCgg — CegCge\^ 
as a function of a in Fig. [31 Note that although 
there are only two points where r = 0, all values of en- 
tanglement, including t = 1 meaning maximal entangle- 
ment, are present in the 'basin of attraction'. Thus we 
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FIG. 3: (color online) The value of the tangle for the states in 
the basin of attraction for different values of a. We notice that 
there are only two points where the tangle is zero, a = ±i. 
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FIG. 4: (color online) The qubit system started in the max- 
imally entangled state (|ee) + \gg))/\/2 and n = 50. (a) the 
entropy of the qubit system, (b) the probability of being in 
the state \gg). (c) the mixed state tangle of the qubit system. 



are observing the time evolution of a generic amount of 
entanglement. To throw further light on the matter we 
show, in Fig. 21 the time evolution of the mixed state 
tangle calculated from p'^ for the maximally entangled 
initial state IV'2) = (|ee) + \gg)) [l3l- Evidently, just 
as the occupation of the initial qubit states collapses and 
revives, so does the entanglement. This phenomenon was 
first noted by Rodrigues et al in a similar context 

Surprisingly, r remains near zero for long periods be- 
tween revivals. Thus, we are dealing with a phenomenon 
which was dubbed the 'death of entanglement' by Yu 
and Eberley [l^ and is the center of much current inter- 
est . Qing et al [2^ have found a similar collapse and 
revival for the same model we have studied here but for 
very different initial conditions. What makes our results 
even more surprising is that the phenomenon occurs for 
a well defined range of initial conditions, namely the 'ba- 
sis of attraction' for all the Nq qubit 'attractor' states, 
and the defining features of these can be generalized to 
an arbitrary number of qubits interacting with the same 
quantum bus. In fact, using the large n expansion of Me- 
unier et al 2l| we have found a 'basis of attraction' for 
all the 'attractor' states in Eq. ([4]) given by 



a if /c is even 



(7) 



if k is odd 



where < |a| < ^J^^^i and the states \Nq,m) are the 
fully symmetrized Nq qubit states, m is the difference 
between the number of qubits in the excited state A^e 
and those in the ground states Ng. 

As noted before, the 'attractor' states are manifestly 



not entangled, but the states in the basin of attraction 
are. Although there is no unique measure of entangle- 
ment for Nq > 2 qubits it is reasonable to assume that 
an arbitrary Nq qubit state in its 'basin of the attraction' 
is gcncrically entangled. Thus the collapse and revival of 
entanglement should be expected to be a generic feature 
of the Nq qubit JCM [13 . 

Interestingly, unlike in the 'two qubit, two cavity' 
model studied by Yonag, Yu and Eberly , in the above 
calculations r decays smoothly to zero with no disconti- 
nuities in the gradient, but does not actually go to zero 
before it revives. That is to say there is no 'sudden death 
of entanglement' 0,[23| and hence there is no need for 're- 
birth' [23]. In fact at t = jtr, when the qubit subsystem 
is in the 'attractor' state, the entanglement is encoded in 
the radiation field. At this time both the qubit-rcsonator 
entanglement and qubit-qubit entanglement vanish. To 
investigate the form this encoding takes we present in Fig. 
[5] the Q function Q{a,t) ~ (a |/0^(i)| a), where p^{t) is 
the reduced density matrix for the radiation field at var- 
ious times. Note that whilst at t = and tr there is 
only one circle which represents a coherent state at 
other times there are two circles. At the interesting time 
t = jtr, when the radiation field is disentangled from the 
qubits, there are two macroscopically different circles on 
opposite sides of phase space so the state of the cavity is 
a superposition of the two coherent states, \a) and |— a). 

Such 'Schrodinger cat' states have been studied by sev- 



eral authors with various perspectives [2i,[2i,[2i. Both 
'Schrodinger cat' states and entangled states may be re- 
garded as particularly 'non-classical' [l^l, while coherent 
states of the field and product states of the qubits are 
regarded as more classical states. As a consequence, this 
fact prompts the following observation: the entanglement 
present in the qubit part of the system at t = is encoded 
in the state of the radiation field, \^{t)) , at t = jtr, which 
becomes highly 'non-classical'. This is demonstrated by 
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shall elaborate on this interesting possibility in a future 
publication [l6j . 

Evidently, whilst we have described the time evolu- 
tion of entanglement in detail only in the 2-qubit limit 
the existence of the 'attractor' states, with a finite basin 
of attraction, for an arbitrary number of qubits implies 
that an oscillatory flow of 'non-classicality' between the 
qubits and the quantum bus, the cavity, is a generic fea- 
ture of dynamics described by the multi-qubit JCM. As 
the properties of this model are relatively readily acces- 
sible, either analytically or numerically, further study of 
the 'collapse and revival' of multipartite entanglement 
dynamics outlined above is clearly called for. In partic- 
ular the effect of decoherence on the persistence of non- 
classicality discovered above remains an open question. 
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edges support from a Royal Society University Research 
Fellowship and the EC QAP project. We thank the ESF 
network AQDJJ for partial support. 



FIG. 5: (color online) Phase space sketches of the Q function 
at six different times when the qubits start in the 'basin of at- 
traction', (a) the time t = 0, where the cavity is in a coherent 
state which is shown by a circle of uncertainty in phase space, 
(b) a time a little after t — 0. (c) the time t — tr/4. (d) just 
before the time t = tr/2. (e) the time 3fr/4. (f) the time 
t = tr, when both the circles have returned to their original 
position. The qubit dipole states are represented as arrows. 
The single arrow at t = corresponds to the spin coherent 
attractor state. 



a comparison of IV'2) hi Eq. depicting the basin of 
attraction, and the analytic result: 



(8) 



For example, if the initial qubit state is not entangled, 
(r = 0), namely a = \^{jtr)) oc \Tct) then the field 
state is a more 'classical' coherent state. However for 
a = where (j) is an arbitrary phase, or a = i?', where 
r is a real number, the qubits are maximally entangled 
(r = 1) and the field is in the 'non-classical' 'Schrodinger 
cat' state |$(|tr-)) « {\a) + |— a)) characterized by a 
Wigner function which takes negative values near the 
origin. In short the 'non-classicality' which was in the 
qubit subsystem at f = is conserved at t = jtr when it 
is in the field subsystem. Remarkably, this also implies a 
new strategy for producing 'Schrodinger cat' states. We 
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